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Abstract 

Multimeson production in pp interactions comprises important background 
for rj, uj and rj mesons production experiments and for the studies of their 
decays planned with WASA detector at COSY. The available information 
about the reactions is summarized and the need for efforts to describe the 
processes is stressed. 
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Direct production of three or more pions in proton-proton interactions, has 
not received proper attention, neither experimentally nor theoretically, de- 
spite the fact that it comprises the main background for rj, r\ and uj produc- 
tion experiments. With a An facility such as WASA, aiming for measurements 
of decays of r\ and rf produced in pp interactions [1] , the understanding of the 
pp —>■ ppmrir reactions becomes very important as they constitute a severe 
background for studies of rj and rf decays into three pions. Those decays pro- 
vide key ingredients for determination of the ratios of light quark masses [2,3], 
since the partial decay widths are proportional to d and u quark mass dif- 
ference squared. In addition precise studies of rj' decays require knowledge 
of pp — > ppinnTiT, pp — > pp7i7i7i7i7i , pp — > ppiTT] and pp —* ppiTTTf] reactions for 
beam energies around r/ production threshold. For more than forty years 
there were only three experimental points available for the cross section 
of pp — > pp7r + 7T~7T° and pp — > pn7r + 7T + 7T~ reactions, all coming from bubble 
chamber experiments [4-6] . Only recently the data base has been extended 
by the measurements of pp — > pp7r + 7r~7r and pp — > pp7r 7r 7r reactions cross 
sections near the threshold by the CELSIUS/WASA collaboration [7]. For 
the remaining reactions there is no data in that energy region. 

The direct production should proceed by an excitation of one or two 
baryon resonances followed by the subsequent decays [8]. For example in 
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the case of three pion production the low energy region a mechanism with 
simultaneous excitation of iV* and A(1232)P 33 resonances is expected to 
dominate. The N* involved has to decay into Nttti and therefore the lowest 
lying Roper (JV(1440)Pn) and iV(1520)Pi 3 resonances could be considered. 

The influence of the resonances can be studied in the invariant mass dis- 
tributions of the subsystems of the outgoing protons and pions. Such studies 
were done for the pp — > ppir + 7T~7T° and pp — > pn7i + 7i + ii^ reactions in bubble 
chamber experiments performed at higher energies (beam kinetic energies of 
4.15 GeV and 9.11 GeV) with up to thousand events [9-11]. However, close 
to the threshold the analysis is not conclusive since the widths of the involved 
resonances are comparable with the available excess energy (Q). Therefore, 
in this case one expects that phase space distribution and the final state 
interaction among the outgoing nucleons provide a reasonable description 
of the observed cross sections. The near threshold cross sections for the 
single meson production via the nucleon-nucleon interaction can indeed be 
quite satisfactory described by such ansatz. For the productions of multiple 
mesons the assumption should hold even for higher excess energies since on 
average the energy available to the pairs of the outgoing particles will be 
lower. That is consistent with CELSIUS/WASA results on pp — > pp7i + n^n° 
and pp — ► pp7r 7r°7r° reactions studied at Q 100 MeV. 

The production mechanism could be studied instead by measuring ra- 
tios of the cross sections for the different charge states. With the lack of 
information the simplest assumption is Fermi model [12], where amplitudes 
for all isospin states are put to be equal. The assumption leads to definite 
predictions for ratios between different charge states of the reactions. For 
example cr(pp — > pp7Y + n^n°) : cr(pp — > pp7r°7r 7r°) : a{pp — > pnn + 7Y + n^) = 
o\ : o"2 : o"3 = 8 : 1 : 10. Resonances in the intermediate state will mod- 
ify the ratios. The effect can be illustrated in the Isobar Model [8] where 
the ratio U\ : <7 2 : o" 3 is 7 : 1 : 25 (assuming AN* intermediate state) or 
5 : 2 : 10 (assuming A^iVJ). Experimentally: (Ji : a 3 at beam kinetic en- 
ergy 2.0 GeV (1:2.53±0.46) [4] and at 2.85 GeV (1:1.59±0.27) [5]. The ratio 
o\ : o"2 = 5.2±0.8 : 1 was measured at lower energy - 1.36 GeV in recent 
CELSIUS/WASA experiment [7]. One trivial modification to the predicted 
ratios close to threshold comes from difference between volumes of the phase 
spaces due to m w + ^ m n o and it amounts to 18% at 1.36 GeV. 

Due to lack of microscopic model calculations, of the same kind as those 
for single meson or for double pion production the semiclassical Isobar Model 
is at present the only option to describe more complicated reactions below 
y/s ~ 5 GeV. The modern version of the isobar model is used as input to 
relativistic ions calculations using transport equations [13, 14]. Reliability 
of the calculations can be tested in simpler cases by comparison with exist- 
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ing calculations for the double pion production [15]. The implementation 
of the resonances in the Isobar Model, their production cross sections and 
decay branching ratios can be evaluated by studying exclusive meson produc- 
tion reactions. Multimeson production in proton-proton interactions provides 
very sensitive test of the parameters. The existing calculations within the 
framework have focused so far on production of dileptons in proton-proton 
interactions with the aim to understand the background for high density nu- 
clear matter probes [16]. The byproduct of such studies were calculations 
of the background for 77 and 77' decays involving dileptons (see for example 
ref. [17]). There is a need to extend the calculations to obtain predictions 
also for multimeson processes. 

2 Background for 77, rj decays 

A feature of 77(77') detection from pp — > ppr)(r}') in the WASA detector is the 
precise tagging by missing mass technique with a resolution of a few MeV/c 2 . 
The resolution of the invariant mass of the decay system is typically consid- 
erably worse. Therefore a figure of merit to describe background from direct 
production process, leading to the identical final state as for an 77(77') decay, 
is given by ps = da b / dp\ M = mjj(rj/) - the differential cross section for the back- 
ground at the 77(77') peak in the pp missing mass p. Figure 1 shows inclusive 
Pb values for 77' derived from the COSY- 11 measurements [18-20]. When 
estimating the background for a given decay channel, quantity p_bA/j (where 
Ap is the resolution in the missing mass) should be compared to a v ^BRi 
(production cross section times branching ratio for the decay mode). The ps 
value depends on the total cross section and on the reaction mechanism. For 
77(77') production close to threshold the background distributions at the edge 
of the phase space are relevant. For both signal and the multimeson produc- 
tion this region of the phase space is strongly influenced by pp final state. 
When approaching the threshold the pb decreases quickly and in addition 
missing mass resolution improves (since it is constrained more by the beam 
momentum resolution) increasing signal to background ratio. The price is 
however a lower production cross section (T v ( v ')- 

Detailed studies of 77 — > nnn decays in pp — * pprj reaction at beam energy 
1.36 GeV by CELSIUS/WAS A collaboration shows that background from 
direct three pion production is 10-20% for the it it + it~ and 5% for the 7r 7r°7r° 
channel in the final data sample [7,22]. This allows for precise study of the 
77 decays providing large number of events is collected. For 77' decays the 
situation is quite different: the three pion decays have branching ratios at 
percent or permil level [23]. In comparison to the 77 meson the 77' production 
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Figure 1: Inclusive differential cross section ps for multipion production 
derived from the COSY-11 data [18-20]. The line is the parametrization 
p B = a(Q/(l MeV))^ fitted to the data: a = 0.64 ± 0.14 [nb/MeV] and 
/3 = 1.66 ±0.08 [21]. 

cross section at similar excess energies is about 30 times lower [18-20, 24- 
26]. Finally the total cross section for multipion reactions increases strongly 
when going from rj to r]' production threshold region, e.g. for pp — > pp-nim 
reaction 50-100 times. In conclusion embarking on the rf decay program in 
pp interactions requires much better understanding of multimeson production 
reactions. 

References 

[1] H. H. Adam et al. [WASA-at-COSY Collab.], arXiv:nucl-ex/0411038. 

[2] D. J. Gross, S. B. Treiman, F. Wilczek, Phys. Rev. D19, 2188 (1979). 

[3] H. Leutwyler, Phys. Lett. B378, 313 (1996). 

[4] E. Pickup, D. K. Robinson, E. O. Salant, Phys. Rev. 125, 2091 (1962). 

[5] E. L. Hart et a/., Phys. Rev. 126, 747 (1962). 

[6] A. M. Eisner et a/., Phys. Rev. 138, B670 (1965). 



4 



A. Kupsc et al 



Multimeson Production 



[7] C. Pauly et al. , Phys. Lett. B649, 122 (2007). 

[8] R. M. Sternheimer, S. J. Lindenbaum, Phys. Rev. 123, 333 (1961). 

[9] G. Alexander et al, Phys. Rev. 154, 1284 (1967). 
[10] A. P. Colleraine, U. Nauenberg, Phys. Rev. 161, 1387 (1967). 
[11] S. P. Almeida et al, Phys. Rev. 174, 1638 (1968). 
[12] E. Fermi, Prog. Theor. Phys. 5, 570 (1950). 
[13] S. Teis et al, Z. Phys. A356, 421 (1997). 
[14] S. A. Bass et al, Prog. Part. Nucl. Phys. 41, 255 (1998). 
[15] L. Alvarez-Ruso, E. Oset, E. Hernandez, Nucl. Phys. A633, 519 (1998). 
[16] A. Faessler et al, J. Phys. G29, 603 (2003). 
[17] A. Faessler et al, Phys. Rev. C69, 037001 (2004). 
[18] P. Moskal et al, Phys. Rev. Lett. 80, 3202 (1998); 
[19] P. Moskal et al, Phys. Lett. B474, 416 (2000); 
[20] A. Khoukaz et al, Eur. Phys. J. A20, 345 (2004). 
[21] M. Zielinski, Diploma thesis, Jagellonian University, in preparation. 
[22] M. Bashkanov et al. , Phys. Rev. C76, 048201 (2007). 
[23] W. M. Yao, et al., J. Phys. G33, 1 (2006). 

[24] A. M. Bergdolt et al, Phys. Rev. D48, R2969 (1993); E. Chiavasa et 
al, Phys. Lett. B322, 270 (1994); H. Calen et al, Phys. Lett. B366, 39 
(1996); H. Calen et al, Phys. Rev. Lett. 79, 2642 (1997); J. Smyrski et 
al, Phys. Lett. B474, 182 (2000). 

[25] F. Hibou et al, Phys. Lett. B438, 41 (1998); 

[26] F. Balestra et al, Phys. Lett. B491, 29 (2000). 



5 



